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1.0 Full second order closure for wall bounded flows. 


1.1 Introduction. 

A full second order closure for wall bounded shear flows is developed, which includes the 
Reynolds stress equations and the equations for all relevant components of the dissipation 
rate tensor. Incompressible and compressible plane flows are considered, but this report is 
only concerned with incompressible flows. 

There are several reasons for the development of a complete second order closure. It can 
be shown that the anisotropy of stress and dissipation rate tensors approaches the same limit 
at the wall, but the derivative of the anisotropy of the dissipation rate is twice the derivative of 
the stress anisotropy at the wall. Another reason is the possibiity of constructing apppropriate 
time scales in the near wall region. The standard second order closure incorporates the 
transport equation for the trace of the dissipation rate tensor and relates the components of 
the tensor via local relations to the trace. The time scale for the destruction of the trace is 
usually modelled using the time scale 



e 

k 


where the modified dissipation rate i is defined by 



with y denoting the coordinate normal to the wall. This is an acceptable model in the region 
close to the wall if and only if the dissipation rate is a nondecreasing function of the distance 
from the wall, because the kinetic energy is of order 0(k ) = y 2 near the wall which implies 
that the second term in the modified dissipation rate is constant. Direct simulations of 
boundary layers and channel flows, however, have shown (see Mansour et ah, 1988) that the 
dissipation rate is a rapidly decreasing function of the wall distance in the viscous sublayer 
and the destruction model using the time scale tq becomes thus a production term in the near 
wall region. This is clearly a violation of realizability for the destruction model. It follows 
that this type of closure model does not represent properly the production of dissipation rate 
near the wall and the sign reversal of the viscous destruction model must make up for this 
deficit in production. The motivation for the modified time scale is the fact that the time 
scale 

_ ^ 

~~ € 

goes to zero as the wall is approached. It will be shown below that the full dissipation rate 
tensor allows a realizable and tensorially invariant construction of a time scale or a time scale 
tensor that reaches a finite and nonzero limit value at the wall. 
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A second point that sets the near wall region apart from the high Reynolds number 
regime of the boundary layer is the growth of the pressure corelations with distance from the 
wall. It will be shown below that the usual split of the pressure correlations into pressure 
transport and pressure rate of strain correlations is not appropriate near the wall, because the 
split correlations grow with different rates and the Taylor series for the original correlation 
involving the fluctuating pressure gradient can be expressed locally up to second order in 
terms of velocity correlations. 


1.2 Exact equations for the dissipation rate tensor. 

Incompressible flows are considered first. Standard manipulations lead to the transport 
equations for the dissipation rate tensor in a Cartesian coordinate system defined by 

e a p = 2vd 1 v' a d 1 v' p (1) 

with expectation denoted by (e Q/ j). The equations can be given in the form 

( d t + (v y )d 1 ){e a p) = <9 7 [i/<9 7 {e Q £) — (u 7 e a ^)] + S^p + S 2 p + 5^ + S*p + S 5 aj} — D a p (2) 


The various source terms are defined as follows. 

= ~{ e P~()^'f( V a) ~ { e ay)dy( v p) (3) 

sJ/.s-(e3j>(a J (»,) + a,K)) (4) 

= -~{ e P'ydy v a) (^) 

sit = ( 6 ) 

Sit = m 

D a p = 4u 2 (d 2 ^v' a d 2 y v'p) ( 8 ) 

Furthermore is the fourth order dissipation tensor defined as 

= 2vd s v' a d y v' fi (9) 


The properties of the source/sink terms on the right hand side need to be established before 
closure expressions can be constructed. 


1.2.1 Equivalent forms of the pressure correlations. 
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The pressure correlations can be given in several equivalent forms. It is instructive to 
split them analogous to the pressure-rate of strain correlations in the stress equations. It 
follows from (6) that 


Sip = -^{d a (d y p'dyv'p) + dp{d y p'd y v' Q )} + B Q p 


where the non-gradient part is defined by 

B a i 3 = ^-{d^p'd^da v'p + dpv' a )) 

The non-gradient part B a p resembles the pressure-rate of strain correlation and shares with 
it the property 

B aa = 0 

It follows that B a p redistributes intensity among the components of the dissipation tensor 
and leaves its trace unaffected. Further splitting of B a p leads to 

B a p = ^-dyid-fp'idaVp + dpv' Q )) - ^-{d^p'idav'p + dpv' Q )) 

The non-diffusive part of B a p contains the Laplacian of the pressure fluctuations which is 
governed by 

= + 2d a v'pdp{v a ) - dlp(v' a v'p) 

This equation has an important consequence: The non-diffusive part of B a p can be repre- 
sented locally in terms of velocity fluctuations. We get 

B a p = —d 1 (d 1 p(d Q v'p + dpv' Q ))+ 

P 

2v{d y v' 6 d 6 vi(d a v' fi + dfiv' a )) + 4v (d y v' s (d a v'p + dpv' Q ))d 6 {v y ) 

We conclude that the non-diffusive and non-gradient part of the pressure correlations does 
not contain a direct influence of the wall and the wall effects can be represented as gradients 
and divergence of a flux. This property very important for the modelling effort. Inspection 
of the local part of B a p shows that it has the structure of the primary production term (5). 
Recasting this part in terms of vorticity and strain rate defined by 

^ or = 


and 


s a p — ^{d a vp -)- dpv a ') 
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leads to 

2u(d y v' s d 6 v' y (d a v l (3 + dj}v' a )) = 4 u(s' afi a' yS s^ s ) - Su{s' Q ^' y u' y ) 
whereas the primary production terms appear as 

Sap = ~^ U (( S a-f S ‘y6 S 6/3) "l" e yotr] ( s 06 S 6'y UJ rj) “b e ^Pr)( s a S S 6y LO r}) ~b e yau e S(3rj { s 6 U u) ) 
The last term can be recast in terms of the Kronecker delta using the tensor relation 

(^7 6 <^ 7/3 ^777 

^ otB ^carj 

&uj 6 ^ijj/3 &u;ri j 


leading to 

$aj3 = ^^ u ii S ay S -f 6 S 6p) “b e yonj {$/36 s 6y UJ rj) “b e 7#rj { S ad 3 6-f U; T ] ) 

+ {’ S 7 /? CJ L a ; 7 } + ( V ^^) “ ~ &<*p{ S 71^7^11)) 

It is apparent that no complete cancellation of triple correlations takes place. The pressure 
correlations appear now as a combination of the divergence of a flux and sources. 

sit = %Ktt + 


where the flux is defined by 

^a/36 “ ”^"( e cj5r7 e u?7<> {&~iP & rj v p) ~b (&yP ^r) V a )) 

or 

Fits = jmp'(bvi + a«0) - - b>(W&,<)) 

The trace of the flux F£g S is not zero but given by 

FL, = ~{djdrf) 

r 

The redistributive sources are given by 

Q% = 2u(d y v' 6 d 6 v' y (d a v' fj + dpv' a )) + 4u(d y v' s (d a v'^ + 0X*)W U 7> 

where the trace of Q v a p is zero. It is noteworthy that one of the components of (d^p'd^v^) 
can be expressed in terms of a component of the dissipation tensor 

4 — {d^p’ div\) = vd 2 (en) + 0 (y ) 

P 
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according to the power series for the pressure given below ( 19 ). For 7 7^ 2 the gradient of the 
wall pressure can be expressed in terms of velocity using the momentum balances. In fact 

~ p d W = vd \y v \ 

holds at the wall y = 0 . 0 . It follows that the wall pressure does not exert a direct influence on 
the dissipation rates if the expansion is carried out to second order. It is clear that modelling 
can be based on the properties of the flux and the redistributive source terms. However, the 
growth rates for the terms in the different formulations of the pressure correlations decide in 
the end their usefulness. This will be investigated in the following chapter. 


1.2.2 Taylor series expansions for the near wall region. 

The near wall region can be analyzed with Taylor series. The coordinate system is 
assumed to be located at the wall and x 2 is the direction of the wall normal pointing into 


the flow field. It is convenient to rename the coordinates and variables as follows: ij = x, 
x 2 = y, X3 = z and v\ = u, v 2 = v, v 3 = w. The velocity components can be expanded with 
respect to the wall normal y 

u(x,y,z,t ) = a 0 +aiy + a 2 y 2 + a 3 y 3 + 0(y 4 ) (10) 

v(x,y,z,t) = b 0 + hy + b 2 y 2 + b 3 y 3 + 0(y i ) (11) 

w(x,y,z,t ) = c 0 + ciy + c 2 y 2 + c 3 y 3 + 0(y 4 ) (12) 

where the coefficients are stochastic functions of x,z y t but not y. They are defined by 


a,(x,z,t) = jf^(O) 

hi*, <) 3 jjajjW 

Cj(x,Z,t)5 jij^(0) 

The noslip condition at the wall implies that 



a 0 = bo = Cq = 0 

( 13 ) 

holds and mass balance 


d a v a - 0 

( 14 ) 

leads to 


dyV 0 = 6l = 0 

( 15 ) 
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and 


~b & z^n — 

for n = 1, 2, ■ * •. The Taylor series for a fixed wall without suction or blowing are therefore 
given by 

u(x,y,z,t ) = a x y + a 2 y 2 + a 3 y 3 + 0(y 4 ) (16) 

v(x,y,z,t) = b 2 y 2 + b 3 y 3 + 0(y 4 ) (17) 

w(x, y , z, t) = ay + c 2 y 2 + c 3 y 3 + 0(y 4 ) (18) 

The expansion for the pressure can be given as follows 

p(x,y, z, t)=p 0 + 2pb 2 y + 3 pb 3 y 2 + p 3 y 3 + 0(y 4 ) (19) 

where p denotes the dynamic viscosity and 

1 d j p 

Pj (x,z,t) = -^(0) 

and the momentum balance normal to the wall was applied to a point at the wall. The 
expansion for the pressure gradient can be shown to be 

(a 2 \ f d x b 2 \ ( pd\ b 3 \ 

d^p = 2p j b 2 1 + 2py I 36 3 1 + 3y 2 j p 3 +0(y 3 ) 

\c 2 J \d 3 b 2 J \pd 3 b 3 J 

which shows that the terms up to first order are proportional to viscosity. The components 
of the Reynolds stress tensor appear in expanded form as 

(u 2 ) = y 2 (o\) + 2y 3 {a I a 2 ) + y 4 (2(a 1 a 3 ) + (a])) + 0(y s ) (20) 

(v 2 ) = y 4 <&!> + 2y*(b 2 b 3 ) + 0(y s ) (21) 

(w 2 ) = y 2 (c 2 ) + 2 y 3 (c lC2 ) + y 1 (2(c 1 c, ) + (c 2 )) + 0(y 3 ) (22) 

{uv) = y 3 (a\b 2 ) + y 4 ((a x b 3 ) + (a 2 b 2 )) + 0(y 5 ) (23) 

The components of the dissipation rate tensor vary near the wall according to 

(«n) = 2i/{(a 2 ) + 4 y(a x a 2 ) 

+y 2 (6(aia 3 ) + 4(a 2 ) + (( d z a x ) 2 ) + ((c?zUi) 2 )) + 0(y 3 )} (24) 

(e 22 ) = 2u{4y 2 (b 2 ) + 12y 3 (6 2 6 3 ) + 0(y 4 )} (25) 

(e 33 ) = ) + 4y{c x c 2 ) 
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+y 2 ( 6 (c x c 3 ) + 4(02) + {(a r ci) 2 ) -f { d z c \ )) + 0 (y 3 )} 


(26) 


(ei 2 ) = 2i/{2y(ai& 2 ) + y 2 (4(a 2 6 2 ) + 3{ai& 3 )) + 0(y 3 )} (27) 

It is apparent from these series that different components grow differently near the wall. It 
follows that the boundary values for the dissipation rate tensor are given by 

(eu)(0) = 2 v{a\) (28) 

(e 22 )(0) = 0.0 (29) 

<e 33 )(0) = 2 u{c\) (30) 

<e 12 )(0) = 0.0 (31) 


It can be shown that the anisotropy of the Reynolds stress tensor is equal to the anisotropy 
of the dissipation rate tensor at the wall and that the normal derivative at the wall of the 
anisotropy tensor of the dissipation rate is twice the normal derivative for the stress ten- 
sor. We consider now the near wall variation of the individual terms in the dissipation rate 
equations. 


1.2. 2.1 Viscous Diffusion. 

The dominant term in viscous diffusion is the normal derivative given by 

dy( u dy( e u)) = 4i^ 2 (6(aia 3 ) + 4(a 2 ) + ((d x a\) 2 ) + {(d z a\) 2 )) + 0(y ) (32) 

d,(vd,{t 22 )) = ViS {b\) + 0(y) (33) 

^(^(ejj)) = 4^(6 (c,c 3 ) + 4(c|> + ((3,c,) 2 ) + {(d.cfj) + 0(y) (34) 

d,(vd,(c 12 )) = 4* 2 (4 (a 2 i> 2 ) + 3 ( 0 , 63 )) + 0(y) (35) 

and they emerge as terms of order unity near the wall for all components. 


1.2. 2. 2 Turbulent diffusion. 

Turbulent diffusion of the dissipation rate component t a p is defined by 

= -a,<qe„„) (36) 
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For the case of boundary layer type flows only the flux normal to the wall is relevant. The 
series expansions lead then to the following expressions for the components of the dissipation 
rate tensor 


dyF l y l — 2j/( 6 2 i (0)) + 0(y 2 ) 

(37) 

d y F 22 =4j/ 3 (6 2 e 22 (0)) + O(y 4 ) 

(38) 

dyFy 3 = 2y{6 2 e 33 (0)) + 0(y 2 ) 

(39) 

d y F 1 y 2 =3y 2 (b 2 d y e 12 (0)) + O(y 3 ) 

(40) 


The expansions show that turbulent diffusion is not of leading order near the wall. 


1.2. 2. 3 Secondary production S^p. 

The interaction of the mean rate of strain and the dissipation rate tensor acts as production 
for the dissipation rates in the same fashion as the Reynolds stress and mean strain rates 
for the stresses. There is however a fundamental difference between this type of production 
for dissipation rates and stresses: It is of leading order for the stress balance but of second 
order for the dissipation rates for high Reynolds number flows. The situation near the wall 
is entirely different. The components of S^p turn out to grow with wall distance as follows 


S, 1 , = -2 y a,(u)(0)a„(e 12 )(0) + 0(y 2 ) 

(41) 

Sh = -Wa 2 „(>,}(0)(&!) + 0(y*) 

(42) 

si j = o.o 

(43) 

Sl 2 = -S^d„{u)(0((bl)+O(y 3 ) 

(44) 


The secondary production terms S^p are not of leading order near the wall, but closed. Hence, 
they need not be neglected. 


1.2. 2. 4 Secondary production S 2 a p. 

The series expansions lead to the following results 

Sn = -y{d„(u)(0)d,(e n )(0) + 2a 2 y (v)(0){e„)(0)) + 0(y 2 ) (45) 

Si, = -2xy , {23,{u)(0)3«<4|) + 8a»,(v)(0 )<*!)} + Otf) (46) 

Si, = -y{d y (u)(0)d,(e„)(0) + 2a^(v)(0)(e M )(0)} + 0(y 2 ) (47) 
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S 2 i 2 = -2uy 2 {2d y {u)(0)((b 2 d xai ) + (a 1 d x b 2 )) + 4d 2 yy (v)(0)(a 1 b 2 )} + 0(y 3 ) (48) 

It is clear from these expansions that the source terms 5^ are of second or higher order for 
boundary layer flows. 

1 . 2 . 2 . 5 Primary production. 

The primary production or vortex stretching terms are the dominant production terms in 


high Re-number flows. For the near wall region they appear in expanded form as 

5 n = + 3 < a iM) + °(y 2 )} (49) 

5*22 = “4i/{y 3 ( 10 ( 62 ) + d x {cii 62 ) + 0(y 4 )} (50) 

Sl 3 = -iv{y{\d,(a,c\) + i(c\b 2 )) + 0(y 2 )} (51) 

S \2 = — 2u{y 2 (d x (b 2 a 2 ) + S(ai& 2 ) + 2(c\b 2 d z a.\) + ( a\Cid z b 2 )) + 0(y 3 )} (52) 


The primary production is not of leading order near the wall, but grows with the same order 
as the secondary production term with wall distance. 

1 . 2 . 2. 6 Viscous destruction. 

Viscous effects can destroy the rate of dissipation and this process is contained in D a p. The 
series analysis shows that the components of D a p are near the wall given by 


Du = 4^{4(a 2 ) + 2((d xai ) 2 ) + 2((d,a 1 ) 2 )} + 0(y) (53) 

£> 22 = 4u 2 {A(bl) + 2iy(b 2 b 3 )} + 0(y 2 ) (54) 

£>33 = 4u 2 {^cl)+2((d x c 1 ) 2 )+2((d z c 1 ) 2 )}+0(y) (55) 

£> 12 = 4i/ 2 {4(a 2 fe 2 ) + y(4(d x aid x b 2 ) + 4(d z a 1 d z b 2 ) + 

12 {a 2 b 3 ) + I2(a 3 b 2 ))} + 0(y 2 ) (56) 

All components of D a p turn out to be of leading order. 

1.2. 2. 7 Pressure correlations. 


11 



The Taylor series for the pressure (19) contains viscous terms which are due to the momentum 
balances. Differentiation of (19) leads to expansions which contain viscous contributions in 
lowest order. The components of the pressure correlations S*p appear as follows 


Sii — 4i^ 2 {6(ai63) + y(2(d x a\d x a,2) + 2(d z a\d z a,2 ) + 4(02^62) + 6(ai<9 x 63)) + 0(y 2 )} ( 57 ) 

522 = 4i^ 2 { 12j/(6 2 63 ) + 0(y 2 )} (58) 

S33 = 4^ 2 {6{cj6 3 ) + y(2(d x Cid x a 2 ) + 2(d z c 1 d z b 2 ) + 4{c 2 d x c 2 ) + 6 (cid x b 3 )) + 0(y 2 )} (59) 


5j 2 — 2i/ 2 {6(0363) +y(2(d x aid x b 2 ) +2{d z a\d z b2) + 12(0263) + 4(62^62) -1 ( a iP3 )) 4- 0(y 2 )} 

pv 

(60) 

The series expansions of the pressure correlations lead to several important conclusion: The 
effect of the pressure correlations in lowest order is local in terms of velocity correlations. No 
Poisson integral appears in lowest order since the wall pressure does not appear in the lowest 
order terms. The split of the pressure correlations obtained in chapter 2.1 leads to a flux 
such that the corresponding source term is strictly redistributive and local in terms of 
velocity fluctuations. The pressure flux is local in velocity in expanded form up to second 
order. Note that there is a cleax advantage for not splitting the pressure correlations in the 
Reynolds stress equations near the wall since the split terms (rate of strain correlation and 
pressure transport) grow with different rates in the viscous sublayer. This is not the case for 
the pressure correlations in the dissipation rate equations. The components of the pressure 
flux emerge for the case of a flat plate boundary layer as 

f , f i 2 = — (ftp's i»;> 

P 

Ay 

F&i = —l(9,p'd ,v' 2 ) + (a jP 'd t vi,)} 

P 

n , 2 = -(vwi 

P 

Ff 22 = — {(ftp'M) - 

P 

which can be analyzed with the aid of Taylor series. We get the following estimates for the 
divergence of this flux near the wall 

9|fT„= 8- 2 <M«ai) + 0(s) 
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= °M 

d,F!,i=&S{hd,c t ) + 0(y) 

^ 2^122 = ~ 4l /2 ( c 2^z c l) + 0(y) 

It follows that 

<9 2 Ffi 2 + d 2 F 3 p 32 = - d 2 (ud 2 (e 22 )) + O(y) 

holds. This proves that the pressure flux terms are of leading order near the wall because 
viscous diffusion is of leading order in this region. 


1.2. 2. 8 Near wall production. 

The production terms proportional to the curvature of the mean velocity profile are respon- 
sible for additional production in the near wall region. The series expansions lead to 

5n=4Ky(^)^{u)(0) + y 2 ((a2)^(u)(0) + 3{a 1 a 2 )^(u}(0) + 2(a 2 6 2 )^(u)(0)) + 0(y 3 )} 

(61) 

S 5 22 = 0(y 3 ) (62) 

*^33 = 0 (63) 

Sfa = vy 2 {{a\)dl yy {u){0) + («iM^»(0))} + 0(y 3 ) (64) 

The near wall production emerges as second order effect near the wall if the boundary layer 
assumptions are satisfied. The component (en) receives all the energy in lowest order. 


1.2. 2. 9 Transport equations in lowest order. 

The series expansions for the source and diffusive terms allow the set up of the transport 
equations for the components of the dissipation rate tensor in lowest order. It turns out that 


all equations are of the same zeroth order. 

<9<(en) = d^{vd^{eu)) — Du + 5^ (65) 

0 = d 7 (i/d 7 (e 22 )) - D 22 (66) 

dt(t 33 ) = d-f( u d~t( e 3z)) — D33 -|- S33 (67) 

0 = d-f(udy(e i 2 )) — D \ 2 + i>i 2 (68) 


These equations are valid near the wall provided none of the surviving correlations vanishes. 
For steady flows the following equations hold then at the wall 

Du — S 31 = dy(udy(eu)) 

D 22 = dy(i/d y (e 22 )) 

D 33 - S 33 = d y (ud 1 (ez 3 )) 

D 12 — S * 2 — d^(vd y (ei 2 )) 

which can be combined with the limit relations for the stress balances to establish constraints 
for the modelled terms at the wall. 


13 


1.3 Closure model for the dissipation rate equations. 


The series expansions for the near wall properties of the dissipation rate equations can be 
used to analyze and to modify closure expressions. In several cases no model expressions 
exist and new models will be developed and analyzed. 


1.3.1 Time scales. 


Several time scales can be constructed for the near wall region with the aid of the dissipation 
rate tensor. First we note that a scale dependent on the wall normal vector can be obtained 
in the form 


k 

n 1 ns{t 1 s) 


( 69 ) 


If n 7 = 8 1 2 it follows that this time scale is given by 


- k 

( £ 22) 

and the series expansions show that both numerator and denominator depend quadratically 
on wall distance. The wall limit is in fact a nonzero value given by 

(dyU 0 dyU0) + ( dyWodyWQ ) 

Hence, a time scale with a nonzero limit at the fixed wall was constructed. This time scale 
avoids the problem associated with the modified (also called homogeneos) dissipation rate 

l = e — 2v{d y \/k) 2 

which may change sign in the flow field. The inverse of another time scale with tensorial 
character using the dissipation rate tensor can be set up as follows 

T ap = ^(e«7)Kvi> _1 (e«/») ( 71 ) 

where denotes the inverse of the Reynolds stress tensor. Conversely is a time scale 

given by 

Ta0 = ~j;( v 'a v> - r )( e ‘t6)~ 1 (v6Vp) ( 72 ) 

where (e Q p)~ l denotes now the inverse of the dissipation rate tensor. 


1.3.2 Turbulent diffusion. 
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The gradient flux model for F° 9 can be given in the form 

Ff=c a ^v' s )d 6 {e aP ) (73) 

where the kinetic energy is denoted by 

*(&*) = \{ v 'a v 'a) 

and the dissipation rate e by 

The constant c 3 has values in the range 0.15 — 0.18. The near wall properties of this model 
follow from the series expansions (10) to (27) as 

dyF" = 0(y 5 ) 


whereas the exact term has a first order variation with respect to the wall distance according 
to (37). Similar discrepancies are observed for the other components. It is clear that model 
expressions developed for the high Re-number regime will not represent the near wall region 
properly. The present model (73) implies that near the wall turbulent diffusion is essentially 
neglected in comparison to the exact term. Inspection of the model (73) shows that there are 
two reasons for its failure near the wall: The time scale approaches zero at the wall and 
the diffusivity is solely determined by the normal stress component {(v 1 ) 2 ) which varies as y 4 
near the wall. The situation can be improved if a composite time scale that approaches the 
scale defined by (69) near the wall is used. The modified closure model is then given by 


F a9 =- 
7 3 


c. 


k 


(v'y s )d S {e<e/3) 


(74) 


The factor 2/3 results from the requirement that the high Re- number limit of the time scale 
must agree with k/e. 

A different model that satisfies all growth estimates and has the correct tensorial and 
dimensional properties can be constructed if the dependence of the turbulent flux on the wall 
parameters n and R ei is taken into account and the near wall model is combined with the 
high .R e -number model. The near wall model is given by 

-{v\e a p)=c'J w2 (R et )[n r) n< : {v'^)]in y (e Q i3) 

It is straightforward to check that it has the same growth rate and the same tensorial prop- 
erties as the exact term. The model represents turbulent transport towards the wall and 
has the form of a convective term. There exist several high Reynolds number models for the 
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turbulent flux of dissipation rate (Hanjalic and Launder, 1972, Lumley, 1978). The present 
model is an analogue of the flux model for the stresses. It is given by 



The combined closure model is the set up as follows 

i k d 

-(w! r 6 a ^)=C , e / u , 2 (ile<)K n cK U <)] Jn 7( C a/3) +C e (l - fw2(Ret))-{ v 'y v 6)-Q^{ e a(3) 


where fwiiRet) denotes a function of the local Reynolds number such that f w 2 goes to unity as 
the wall is approached and to zero in the turbulent zone. Furthermore, The function f W 2 {Ret) 
should be nonnegative and it should not modify the dependence on the wall distance for the 
near wall model. It follows that f w2 must be an exponential function of the Reynolds number 
given by 

fw2(Ret) = exp[-(^-) 2 } 

U e 2 

where R e2 is a constant measuring the range of influence for the near wall model. This 
function of the Reynolds number has the well known property that all its derivatives at zero 
vanish. Hence, it does not modify the growth rate of the near wall model. 


1.3.3 Secondary production S 2 a p. 

The properties of Sip in the high and low Re- number limits will be considered first. Local 
isotropy requires that 

holds. It follows that the secondary production terms Sip are for high Re-numbers given by 

lim Sip = -^6 a ped y (v y ) (75) 

ite— +oo y 

The divergence of the mean velocity is zero for incompressible flows and it follows that the 
high Re-number form of the secondary production terms can be neglected. The near wall 
variation of the secondary production term 5^ according to (45) can be regarded as 

s u = -(u)d x (e u ) - 2(e n )d y (v) + 0(y 2 ) 

and it follows that this term is at best of the same order of magnitude as the mean convection 
term. Similar relations hold for the other components and we conclude that the secondary 
production terms Sip can be neglected. 


1.3.4 Primary production and viscous destruction. 
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The closure model for these terms is fundamentally different in the high and the low Reynolds 
number limits. Both limits need to be considered and the corresponding closure models must 
be merged to cover the range of Reynolds numbers from zero to infinity. 

High Reynolds number limit. 

The order of magnitude estimates for primary production and viscous destruction at high 
Re-numbers shows that they are of leading order, but their difference is of the same order 
as the secondary production term Sip. It follows that they should not be treated separately 
but their difference should be modelled as function of the available information. The model 
consistent with second order moments is in general given by 

Sip - D Q p=-V a p({e Su; ), {v'svl), ds{v u,)) (76) 

The dimensionless and symmetric tensor ^ a p should represent both production due to the 
interaction of vorticity and strain rates and the destruction due to viscous effects. If we 
impose the condition of local isotropy on this model we get the following variant 

d y (vs)) - ¥$((«« w ), d 6 (v w )) 

where the first part represents the productive and the second the destructive contribution to 
the model for the difference of and D a p. The model can be set up to be consistent with 
the standard expressions for the trace equation (see Launder, Reece and Rodi, 1976) 


^ a{3 — 




fa) -C e 2 


€ 


and the time scale is given by 

k 

e 

The closure model 


Sip - D a p= - c t i^ a pj{v\v' 6 )d^{v 6 ) - c (2 (e a p ) | (77) 

emerges. It is, however, not applicable to wail bounded flows since the destructive part of 
this model becomes singular as the wail is approached. This deficiency can be corrected 
either by defining a time scale that does not go to zero at the wall or by merging low and 
high Reynolds number closures with a Reynolds number dependent function such that the 
singularity is removed. Finally, we note that this model is not necessarily positive definite 
because the productive part is not positive definite, a property shared with the exact 
term. Refined closure models can be constructed using tensorial time scales introduced in 
chapter 1.3.1. 

Low Reynolds number limit. 
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The dependence on the wall distance is the deciding property as the wall is approached. We 
recall that the primary production terms decay with wall distance, according to chapter 2.2.5, 
as follows 

0(S 3 U ) = 0(Sj 3 ) = y 
0(Sf 2 ) = y 3 

0(S 3 12 ) = y 2 

whereas the viscous destruction terms are all of leading order 


O(D Q0 ) = 1 


according to chapter 2.2.6. It follows that they must be modelled separately in the near wall 
region. The first step in the construction of the low Reynolds number version is the analysis 
of the near wall properties of the high Re-number model. It follows from (77) that 

0(^* (1) ) = y 

which is at variance with the detailed decay laws for 5 3 ^ given above for a = (3 = 2 and 
a = 1,/? = 2. However, since the primary production is not of leading order near the wall 
it would be acceptable. A more serious problem arises in the destructive part If we 

construct a time scale such as (69) for the near wall region we avoid the singularity, but it is 
not possible to satisfy the decay laws for the destructive terms D a p. It follows that the time 
scale in any closure of the form given by cannot be a scalar but has to be a tensor of 
rank two (or higher) with positive eigenvalues. It is not difficult to construct a closure model 
such that the decay law is satisfied. For instance, the model 


D a p= - c (2 r 6 


dxsdx v dxfidxrj 


with (69) as time scale possesses the correct tensor properties and the correct decay law as 
the wall is approached. However, it is unacceptable as closure model because it produces 
in regions where the second derivatives are all positive a second order pde with negative 
diffusivity. The initial/boundary value problem for such an equation is not well posed and 
the numerical solution futile. Furthermore, this model would not be realizable, because the 
limit (e a p) — ► 0 does not imply that the second derivatives of (e a p) go to zero. It follows that 
the closure model for the viscous destruction of the dissipation rate components must be of 
the form 

Da/3 = F'a0(( € y r )), ( v 7 v n )j ^c( C 7»f)’ ^c( V y V i))' Ret) 
where R et denotes the turbulent Reynolds number defined by 



tv 


(78) 
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The example cited above for the destruction of dissipation rate showed that derivatives lower 
that second should be used. Several closure models can be set up that are tensorially consis- 
tent, have the correct decay law near the wall and avoid the stability problems incurred by 
negative diffusivities. We note two of them, first 


Do,0 — c e2 


(jod) 

T 




d( e ar?) 0(6/0 
dx y dxt 


(79) 


and 

n — — r P-L — r i/er> n — \ /om 

D a0 - c e2 ^ Ct 3 ven y n 6 g )^( g ) (SO) 

as second variant. Note that n Q denotes the unit normal vector of the wall pointing into 
the flow field. The time scale r is given by (69) to avoid a singularity at the wall. The two 
models look very similar, but inspection of the time scale (69) shows at once that the first 
model is most likely unstable. Suppose that kinetic energy k and dissipation rate component 
n y ns{€ y s) = (^ 22 ) are in equilibrium, then assume that this equilibrium is disturbed, say the 
kinetic energy is reduced by a small amount. It follows from the fact that the second part in 
the first model is quadratically proportional to the kinetic energy that the rate of destruction 
of the dissipation rate components is decreased by the disturbance and consequently are the 
normal stresses further decreased and the equilibrium state is not recovered. Hence is this 
model unstable. It follows that the second model is the preferred one. 


Merger of low and high Re-number models. 

The closure models for the destruction term (79) or (80) valid for the limit R e — > 0 and 
the destructive part of (77) valid for the high Re number limit must be merged together to 
produce a model valid throughout the flow field. Suppose both limit expressions have decay 
laws near the wall that do not need change via the function weighting them according to the 
local Reynolds number (or any other function propertional to the distance to the wall). The 
we need a weight function that does not change the dependence on wall distance near the 
wall. This implies that we must find a function which has zero derivatives at zero. It is well 
known that the exponential function 


f(y) - ex P(~^) 

has the desired property, in fact, it is infinitely often differentiable but not analytic at zero 
(its Taylor series is identically zero at y — 0). Hence, we can establish a low Reynolds number 
function 

/(Re) = ex P[“(^) 2 ] 

which is zero at zero Re-number and unity at infinite Re-number. The constant R° determines 
the range of Reynolds numbers for which the function is close to zero. Other functions 
have been proposed that are proportional to some power of the wall distance and change 
therefore the decay law. Several models have been proposed for the functions f \ , multiplying 
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the production of dissipation rate, and / 2 , multiplying the destruction of dissipation rate. 
Vandromme et al. (1983) suggested 

fi = i-o 

and 

Up 

/ 2 = 1.0-0.22exp{-(^) 2 } 

based on the original model of Hanjalic and Launder (1976) and obtained good agreement 
with measurements in flat plate boundary layers. Recent developments surveyed by Launder 
(1989) use the second and third invariants of the Reynolds stress tensor to represent the wall 
influence on production and destruction processes. 


1.3.5 Pressure correlations. 


It was shown in chapter 1.2 that the pressure correlations can be split into transport and 
source terms such that the source terms are strictly redistributive and have no effect on 
the trace of the dissipation rate tensor in analogy to the pressure-strain correlations for the 
Reynolds stress tensor. It follows that they must redistribute intensity among the components 
of the dissipation tensor. Kolmogorov’s hypothesis of local isotropy requires that the dissi- 
pation rate tensor approaches its isotropic form as the Reynolds number approaches infinity. 
It follows that a return to isotropy model given by 

- 7«M - (W) 

would satisfy this condition. The open question is the time scale. The scale 

_ k 

e 

is the obvious choice for the high Re-number limit, but the model becomes incorrect as the 
wall is approached because r -1 goes to infinity with y~ 2 . Modification of the time scale 
according to (69) solves this problem and 


<%„= - c, - f (SO) 

emerges as nonsingular variant. The transport part of the pressure correlation can be mod- 
elled in terms of the viscous diffusion terms because the wall limits indicate this form. The 
model 

where the tangential unit vector is defined by 


_ Mm 


( 82 ) 
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and the binormal by 


satisfies the near wall properties for the components F[ 12 and F ^ 32 and neglects the compo- 
nents F 222 an d F[. 22 - The effect of the diffusive part of the pressure correlations is therefore 
inhibition of the viscous diffusion near the wall for the diagonal components corresponding 
to motion parallel to the wall. The effect on shear component and the diagonal component 
corresponding to the motion normal to the wall are neglected. 


1.3.6 Near wall production. 

The near wall production requires in general flows a closure model. However, for the near 
wall region in boundary layers expressions can be given that are exact in the wall limit. 

Sn= - 4vd 2 (v' 1 v' 2 )dl 2 {v l ) (83) 


*$22 — ^33 — 0.0 (84) 

S\ 2 = - ud 2 {v 2 )dl 2 (v{) (85) 

The model assumption is essentailly the assumed validity of these expressions for finite dis- 
tance from the wall. Since all expressions axe proportional to the laminar viscosity, quick 
decay with increasing wall distance can be expected. 
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2.0 Reynolds stress and complete second order models. 


The usual closure models based on the trace of the dissipation rate tensor will be considered 
first. Several versions of the Reynolds stress model are available and the most relaible stress 
model (Vandromme et al., 1983) will be evaluated. It will serve as test bed for the complete 
second order closure. 


2.1 Stress equations. 

The balance for the Reynolds stress components is given by 

= -M «”>;>¥ + + £» 

_ (p^v'pv'J - 8 ay {p'v'p) - 8p y {p' v' a )) - (p)e aj3 


dx. 


( 86 ) 


where the dissipation rate tensor is defined by (1). The split of the pressure correlation into 
pressure transport and pressure rate of strain correlations can be shown to be inappropriate 
near the wall. It follows from (19) that the Taylor series for the pressure gradient is given in 
terms of velocity derivatives at the wall up to second order. Hence, it is possible to represent 
the Taylor series for correlations involving the pressure gradient in terms of local velocity 
correlations at the wall up to second order. If the pressure correlations are split as in (86) 
the pressure fluctuation itself appears and the solution of the Poisson equation introduces 
the well known integral contributions. 

The standard second order closures employ the equation for the trace of the dissipation 
rate tensor 

es-(e aa ) 

which follows at once from (2) 

(df + (v y )dy)e = d y [ud y e - (v' y e) + F?] + S 1 + S 2 + S 3 + S 5 - D (87) 

The source terms are defined as follows. There are two groups of secondary production terms 

S 1 = ~(e ay )d y (v a ) (88) 

and 


S 2 = -(udsv' a d y v' a )(d y (vs) + d 6 (v y )) 
The primary production term is given by 

S 3 = -(e„, d y v' a ) 


(89) 


(90) 
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and the pressure correlations can be given in the form 


S 4 = (d^p'd^) (91) 

It was shown in chapter 1.2.1 that the pressure correlations can be represented as the sum of 
the divergence of a flux and a redistributive source which has zero trace. Hence we get 

S 4 = d 8 Ff (92) 


where the flux is defined by 


d v v' Q ) (93) 

or 

F! = -jidyp'd^l) (94) 

The near wall production of dissipation rate is proportional to the curvature of the mean 
velocity profile 

S 5 = - 2iy(v' s d y v' a )d^ 6 {v a ) (95) 

and the viscous destruction of dissipation rate is given by 

The properties of the source terms have been established in chapter 1. and we observe that 
none of the source terms is closed in contrast to the equations for the dissipation rate tensor. 


2.2 Standard second order closure model. 


The development of the full second order closure was based in a systematic way on 
existing closure schemes. The standard second order closure using the trace of the dissipation 
rate tensor set up by Vandromme et al. (1983) was an important stepping stone and can be 
applied to test closure models for the various processes governing the dissipation rate tensor. 
The model is given by 




+ <»K> 


d, r 7 


,g. , d , dWp) 

+ 3’a/3 + -X {p- 


dx. 


dx. 




( 97 ) 


s y V/O/ 

The pressure rate of strain correlations are modelled by 

*.«= - 9/,(s,)Wj((»>;) - !*„,*) - /,(R.)W{ 2 jp(j’^ - \ksP) 
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( 93 ) 


— n - (fl “'’ - 3 S ^ P) - -^r k( -d-x 7 + + *“' } 

as the sum of return to isotropy, fast response and wall contributions, where 


r> / / lxd(v0) / / 


(99) 


and P = l/2P aa and 


-i«i“ - (.»« 


The wall contribution is given by 


5{v q ) 


♦JiWW - 3^o») + 4(^ - D./>) + <**(-§^ 


( 100 ) 


3W _ 2 ( 


3:C„ 3~ af ' dx. 


■)}■ 


e y 


(101) 


The dissipation rate tensor must be modelled in terms of its trace and the anisotropy of the 
stress tensor 

(e Q0 )=j(fs(v' a v' fi ) + (1 - f s )\^k) (102) 

where f 3 depends on the local Reynolds number. The turbulent flux is modelled by 


F<*0 7- - c,t(p) e { v 'i v 'f>) Q XB 

(104) 

The constants are given by c st = 0.25, ci = 1.5, C 2 = 0.4, c' x 
c 5 = 0.0041. The low Reynolds number functions are given by 

= 0.1597, 4 = 0.0133, 

/.(*.) =«p( 1 + b</30 ) 

(104) 

and 


fp(R e ) = tanh(R e /5 0) 

(105) 


The equation for the trace of the dissiaption rate tensor requires closure for all relevant 
processes on the right hand side. Following Vandromme et al. (1983) the model 




-Cfi{p)h{Rt)j\e - 2 1 + 


k dx-fdxs dx-fdxs 


dx Q 

(106) 
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The low Reynolds number functions are given by 


/„ = [!- exp( — /„i i? y )] 2 (l + 



( 107 ) 


and 


where 


fl 1 + 8 (/„ + 10- 10 ) 3 

h = 1-0.22 exp(-S) 



(108) 

(109) 

(110) 


and /„i = 0.0165, f„ 2 = 20.5. The performance of this closure model was tested in fiat plate 
boundary layers. The results are contained in fig.l to fig. 8. The mean velocity in fig.l is in 
good agreement with the law of the wall. The turbulent Reynolds number and the damping 
function f s (R e ) in fig. 2 and fig. 3 prove that the boundary layer is fully developed. The 
Reynolds stress components in fig. 4 to fig. 8 show the expected distributions with maximal 
values in reasonable agreement with the experiments. 
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Figure 1 Variation of normalized mean velocity 
with wall distance in wall units 
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Figure 3 Damping function (F s )for dissipation rate 
term in Reynolds stress equation vs. y + 
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Figure 7a Variation of normalized sheer stress u’v' with 
wall distance in boundary layer thickness 
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Figure 7£Variation of normalized sheer stress u v with 
wall distance in wall units 
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Figure 8 Variation of normalized turbulent kinetic energy 
with wall distance in wall units. 
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2.3 Full second order closure model. 


The development of the full second order closure was prepared in the first chapter which 
contains the properties of all transport equations for second order moments derived from 
series expansions with respect to wall distance. The results obtained with series expansions 
are called growth laws since they describe in first or higher order the growth of correlations 
with wail distance. The challenge is now to construct a closure model that satisfies all tensorial 
and realizability conditions and the growth laws. 


2.3.1 Closure model for the stress equations. 

The closure model for the stress equations follows closely the model developed by Van- 
dromme et al. (1983). The only difference is that the model (102) for the components of the 
dissipation rate tensor is not used. The stress equations appear then according to chapter 
2.2 as 


-Clfp{Re){p)^((v' a v'p) - \s a pk) - f p {Re){p){ ^J^ (P a (3 ~ \ 8 apP) ~ ~ ( D <*0 ~ \ 8 ^P P ) 
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30c 2 -2 ,,d{vp) d(v a ) 


55 


■*( 


dx f 




, a d(v' a v' 0 ) 
*dx y ifi dx 1 




( 111 ) 


where p a 0 and D a p are defined in (99) and (100) respectively and the wall contribution 
to the pressure correlations is given by (101). The low Reynolds number functions are all 
established in chapter 2.2. 


2.3.2 Closure model for the dissipation rate equations. 


The transport equation for the dissipation rate tensor 

+ ( v f)~Q^)( e °‘P) = ~^ L + Sip + + Sip + S*p + Sip — D a p (112) 

requires closure for the turbulent flux contained in the total flux F a p~, and the source terms 
Sip, Sip, S^p, S^p and D a p. 
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Model for the turbulent flux: This model follows the suggestion of Vandromme et ah 
(1983) for the turbulent flux of the trace of the dissipation rate tensor 


•F a/)~) — ^ 


dfcap) 

dXy 


+ C c{p) 


k 





d(ta£) 

dxs 


(113) 


The only modification is the use of the time scale (69) to improve the behaviour near the 
wall. The value for the constant c £ =0.1 is consistent with the six-equation model. A more 
sophisticated model was developed in chapter 1.3.2. 

Model for the source S 2 . j : The secondary production can be neglected near the wall 
according to chapter 1.2. 2. 4. 

Model for the source S The primary production is not of leading order near the wall 
but varies like the closed production term in the near wall region according to ch.1.2.2.5. 
The present model utilizes this property and the analogy to the stress transport exploited 
in the equation for the trace (106), where the high Re-number part of the model for the 
difference between the leading terms having the character of a production term is modelled 
proportional to the production of kinetic ernergy 


S lp = - C f 3 (p)i( € 0y^~ + ( 114 ) 


The constants are given by c £ i = 1.45 and c £ 3 = 1.0. 

Model for the source term S*pi The pressure correlations were shown to be of leading 
order near the wall for all dissipation rate components except €22 ( see ch.1.2.2.9). They can 
be split into diffusive and reditributive source terms (ch.1.2.1). The diffusive part is assumed 
to be represented by the closure for the turbulent flux (113). The model for the redistributive 
source is analoguous to the return to isotropy model for the stress transport equations. It is 
given by 


~ c «4 (p) , (( e a/j) — <>^otj9 e ) 


(115) 


with c {4 = 12.5. 

Model for the destructive term D a p: It was shown in chapter 1.3.4 that the high Re- 
number model for the difference of primary production and viscous destruction becomes 
singular as the wall is approached. Furthermore, the growth law for the viscous destruction 
term D a p implies that a model similar to the high Re-number case 


Dafl Ct2 (p) 


(fog) 

T 


with a scalar time scale is impossible. It follows that a tensorial time scale m ust be con- 
structed to conform with the growth law. The present closure is a composite expression 
containing the high and the low Re-number models 


-Da£=c £ 2i(p)(l 




2 

nrjTlJ, 


djtory) djtpj) 
dx,j dx $ 
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+Ce 22 (p)n a n l} n y ns( 1 - f 3 (R e t))j(e lS ) + c f 23 (e Q /?) (116) 

where the preliminary values for the constants are c e2 i = 0.32 • 10 — 4 , c t2 2 = 26.25, and 
c ( 23 = 12.5. The low Re-number functions are set up as follows: 

20 5 

fu(R et ) = (1 - exp[-0.0165R y ]) 2 (l - — -) (117) 

Met 

where 

Ry 

denotes the dimensionless wall distance and 

fs(Ret) = tanh(0.004R et ) 

The turbulent Reynolds number is defined by 

Ret = ~ 

ev 


(113) 



2.3.3 Preliminary results for the complete second order closure. 

The system of nine parabolic differential equations was tested in reduced form by pre- 
scribing the profiles for mean velocity and the Reynolds stress components which were ob- 
tained with the six equation model discussed in chapter 2.2. The numerical solution for the 
remaining equations for the disssipation rate tensor was carried out and convergence was 
achieved after a few hundred steps. The results axe presented in fig. 9 to fig. 16. The figures 
contain also as broken line the dissipation rate components deduced from the local relation 
suggested by Launder and Reynolds (1983) and modified by Lai and So (1990) 


where 


+fw(Ret) 


( e or/?) — g(l fw{Ret)^al3^ 

t ( V a V p) + ( V a V '-,) n -t n 0 + ( u ^)n 7 n a + U a Tl pH y n S {v^Vg) 

k 1 + ^t n 7 n «K u «> 


fw(Ret) - eX P[“(^) 2 ] 


(119) 


The dissipation rates in the figures are normalized with the wall variables v and u T 

{taf})'' 


= 
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with u r = y/r^Jp denoting the wall shear velocity. The symbols in the figures represent the 
results of the numerical experiments carried out by Mansour, Kim and Moin (1988). 

The component in fig. 9 (logarithmic scale) and fig. 10 (linear scale) shows a negative 

gradient near the wall which implies that the time scale 


T 


-1 





2 


would change sign and the destructive term would become a production term violating real- 
izability. The prediction of in the outer (fully turbulent) part of the flow field is too small 
and the near wall part appears too large compared to the direct simulation, but the profile 
shape is in good agreement with the numerical experiment. The component €22 fig- 11 and 

fig. 12 reasonable agreement between the full second order closure and the direct simulation, 
but the local relation of Launder and Reynolds overpredicts the component by a factor of 
three. The prediction of the component ef 3 in fig. 13 and fig-14 shows a similar behaviour as 
the component e in fig. 9/10. The shear dissipation ef 2 * n fig-15/16 shows overprediction by 
the full second order closure and underprediction by the local relation near the wall whereas 
the outer part is in good agreement with the direct simulation data. 


2.3.4 Conclusions. 

The results presented lead to several conclusions. It is clear from the theoretical develop- 
ment that only the full second order closure offers the tools to construct the appropriate time 
scales in the near wall region of a turbulent boundary layer. The growth laws for the various 
correlations appearing in the stress and dissipation rate balances limit severely the model 
expressions and indicate that composite models for the high and the low Reynolds number 
regimes must be established. The model discussed in this chapter produces good results if 
velocity and stress components axe held fixed and this indicates that the model expressions 
are consistent with the direct simulations. However, the stability of the closure needs to be 
investigated and this part may lead to modifications of the present version of the full second 
order closure. This part of the project is currently under way. 
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Figure 9 .-Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. A Launder & Reynolds. 
— :Present Work. 
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Figure 10 :Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. ▲ Launder & Reynolds. 
— :Present Work. 
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Figure 11 :Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. ▲ Launder & Reynolds. 
— :Present Work. 



Figure 12 : Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. ▲ Launder & Reynolds. 
— :Present Work. 
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Figure 13 : Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. A Launder & Reynolds. 
— :Present Work. 



Figure 14 .Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. ▲ Launder & Reynolds. 
— :Present Work. 
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Figure 15 Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. A Launder & Reynolds. 
— .Present Work. 



Figure 16 : Variation of Psuedo 9 equation model of normalized 
dissipation rated with wall distance. O : Results from Mansour, 
Kim, and Moin direct simulation results. A Launder & Reynolds. 
— :Present Work. 


37 



References. 


Hanjalic, K. and Launder, B.E. (1976), ’’Contribution towards a Reynolds stress closure for 
low-Reynolds number turbulence”, JFM 74, 593. 

Lai, Y. G. and So R. M. C. (1990), ”On near-wall turbulent flow modelling”, J. Fluid Mech. 
221 . 641 

Launder, B.E. and Reynolds, W.C. (1983), ’’Asymptotic near-wall stress dissipation rates in 
a turbulent flow”, Phys. Fluids, 2§, 1157 

Launder, B.E. (1989), ’’Phenomenological modelling: Present and future”, Whither Turbu- 
lence? or Turbulence at the Crossroads Conf., Cornell Univ. 

Lumley, J. L. (1978), ’’Computational Modeling of Turbulent Flows”, Adv. Appl. Mech. 18 . 
123 

Mansour, N. N., Kim, J. and Moin, P. (1988), ’’Reynolds- stress and dissipation-rate budgets 
in a turbulent channel flow”, J. Fluid Mech. 192 . 15 

Shima, N. (1988), ”A Reynolds stress Model for Near- Wall and Low-Reynolds-number Re- 
gions”, Trans. ASME J. Fluids Engng. 110 . 38 

Vandromme, D.D., Ha minh, H., Viegas, J.R., Rubesin, M.W. and Kollmann, W. (1983), 
’’Second order closure for the calculation of compressible wall bounded flows with an implicit 
Navier-Stokes solver”, 4th Symp. on Turbulent Shear Flows, Karlsruhe, Germany. 


38 


